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We investigate the interplay of collective and chaotic mo- 
tion in a classical self-bound iV-body system with two-body 
interactions. This system displays a hierarchy of three well 
separated time scales that govern the onset of chaos, damp- 
ing of collective motion and equilibration. Comparison with a 
mean-field problem shows that damping is mainly due to de- 
phasing. The Lyapunov exponent, damping and equilibration 
rates depend mildly on the system size N. 

PACS numbers: 24.60. Lz, 24.30. Cz, 05.45. Jn 



Self-bound many-body systems like the atomic nucleus 
display regular, collective phenomena as well as chaotic 
behavior [Qj2|. The giant dipole resonance, for example, 
constitutes a regular oscillation of the protons against the 
neutrons albeit it is excited at energies where the spec- 
trum displays fluctuations which are typical for chaotic 
systems ||. This interesting interplay of collective and 
chaotic motion and the effects of chaotic dynamics on 
the damping and dissipation of nuclear excitations is a 
matter of intense research Q Q but still not fully un- 
derstood. 

Several authors P,p|,p|, ^0p^ ] have addressed equilibra- 
tion and damping of collective motion by coupling a 
slowly moving collective degree of freedom (typically the 
wall of a container) to fast moving independent parti- 
cle degrees of freedom. Within these models, collective 
degrees of freedom are neither constructed from single- 
particle coordinates nor the result of a self-consistent 
mean field. In the case of driven walls this may lead to 
energy conservation problems. Bauer et al. [H treated 
a more realistic model where the collective motion was 
coupled self-consistently to the single-particle motion and 
found the coexistence of undamped collective motion and 
chaotic single-particle dynamics. This observation is in- 
teresting with view on low-lying collective excitations 
but cannot apply to the case of the giant dipole reso- 
nance. Furthermore, the classical models discussed in 
the literature suffer from the absence of a two-body force 
and rotational symmetry. In this work we are particu- 
larly interested in the many-body aspects of the problem 
and the role of the two-body interaction opposed to the 
mean-field. To this purpose we study a classical self- 
bound many-body system with the following character- 
istics: two-body surface interaction; chaotic dynamics; 
collective mode constructed from single-particle degrees 
of freedom. 

Let us consider the model Hamiltonian 
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where fi is a two-dimensional position vector of the i-th 
nucleon and pi is its conjugate momentum. The interac- 
tion is given by 



V{r) = 



for r < a, 
oo for r > a. 



(2) 



The particles thus move freely within a convex billiard in 
2iV-dimensional configuration space and undergo specu- 
lar reflection at the wall. This corresponds to the basic 
picture of nucleon motion in a nucleus: they move in 
a common flat-bottom potential, and interact mainly at 
the surface. The Hamiltonian achieves these fea- 

tures while retaining the two-body interaction. Total 
momentum P, energy E and angular momentum L are 
conserved. Classical phase space structure is indepen- 
dent of energy since a scaling pi — > apl simply leads to 



a rescaling of E 



2 E, P -» aP, and L -> ah. We 



set m = a = H = 1 . Thus, energies, momenta and times 
are measured in units of H 2 /ma 2 , h/a and ma 2 /h~, re- 
spectively []. In what follows we fix the total momentum 
P = 0, the angular momentum L — and set the energy 
E = N. 

The N = 2 system is integrable and chaotic dynamics 
may exist for N > 3. To study the classical dynamics of 
the ./V-body system we compute the Lyapunov exponent 
using the tangent map []l5| . The time evolution of the 
TV-body system is done efficiently by organizing collision 
times in a partially ordered tree and requires an effort 
0(N \n N) We draw several initial conditions for 

fixed N at random from the (E — N, P — 0, L — 0)-shell 
and follow their time evolution for several hundred Lya- 
punov times. This ensures good convergence of the Lya- 
punov exponents. The results of several runs are listed in 
Table Q (mean values and RMS deviations). For N = 3 
we traced 7 x 10 4 trajectories, all of them being unstable 
with respect to small initial deviations. Thus, we ex- 
pect more than 99% of phase space to be chaotic. For 
larger N we traced fewer trajectories and have less statis- 
tics. However, all followed trajectories have positive Lya- 
punov exponents. This suggests that the dynamics of the 
./V-body system is dominantly chaotic. We note that the 
iV-body system possesses marginally stable orbits corre- 
sponding to configurations where N — 2 particles are at 



1 We introduced Ti to set an energy scale. 
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rest. However, these configurations are of zero measure 
in phase space. The reliability of the numerical evolu- 
tion was checked by comparing forward with backward 
evolution. Moreover, total energy, total momentum and 
angular momentum were conserved to high accuracy. As 
a further check we computed the Lyapunov exponents us- 
ing an alternative method |l7j and found good agreement 
with the results displayed in Tabic jjj. 

We next consider the time evolution of collective mo- 
tion. To this purpose we define a set of initial condi- 
tions (i.e. phase space points) of the iV-body system 
that lead to collective motion corresponding to the gi- 
ant dipole resonance. In passing we mention that one 
could also use collective coordinates introduced by Zick- 
cndraht fl§f| or consider motion close to invariant man- 
ifolds of the rotationally symmetric many-body system 
fl9f . Let us draw the momenta pi at random from uni- 
form probability distributions that vanish outside of the 
domains (p x ,p y ) G [-Aq x , Aq x ] x [±q - Aq y , ±q + Aq y ] 
, with + or — sign for particles i = l,...,N/2 (e.g. 
protons) and particles i — N/2 + 1,...,N (e.g. neu- 
trons), respectively. We rescale the energy to E = N. 
The initial positions fi , i = 1, . . . , N — 1 are drawn at 
random to lie inside a circle of radius a/2; the posi- 
tion of the Nth particle is chosen such that the total 
angular momentum vanishes; the center of mass is sub- 
tracted. For Aq x w Aq y <C q ~ ^j2mE/N one obtains 
initial conditions that correspond to the motion of all 
protons against all neutrons and therefore may be iden- 
tified with the dipole giant resonance^. In what follows 
we use Aq x — Aq y w g/10. 

We are interested the time evolution of the dipole- 
moment 

2 AT/2 
i=l 

Fig. [I] shows the time dependence of the dipole-moment 
for a system of size N = 220 and one collective initial 
condition. The longitudinal component D y exhibits a 
few damped oscillations that eventually turn into erratic 
fluctuations of small amplitude. On the same time scale 
the transverse component D x increases in amplitude from 
approximately zero and turns into similar erratic fluctua- 
tions. These fluctuations decrease with increasing system 
size N , indicating that they are of statistical nature. The 
first two periods of the time evolution can be fitted to an 
exponentially damped harmonic oscillation yielding the 
period r and the damping rate 7. The results of sev- 
eral runs (mean values and RMS deviations are listed in 
Table [n|. 



2 Our choice of zero angular momentum is justified in the 
framework of classical mechanics. 



For the equilibration we study the distribution of 
single-particle energies ej , j = 1 . . . N during the time 
evolution of the dipole mode. The moments I v = 
N^ 1 €j of this distribution may be compared with 

the corresponding equilibrium values obtained from a 
Boltzmann distribution p {B) (e) = l/[fcTexp (e/kT)] at a 
temperature kT = E/N = 1, i.e. I^ B) = v\{kTf . Fig. | 
shows the first nontrivial moment I2 averaged over many 
runs to reduce statistical fluctuations. This moment ap- 
proaches its equilibrium value exponentially fast at a rate 
a (listed in Table ||) . The saturation at long times seems 
to be due to the finite number of particles. Once the 
system has equilibrated, the single-particle momenta are 
Maxwell-Boltzmann distributed (See Fig. ||). We also 
computed the Lyapunov exponent A of the many-body 
trajectory corresponding to the collective dipole motion 
and listed the results (mean values and RMS deviations) 
of several runs in Table ||. The Lyapunov exponents 
agree well with the previously computed ones listed in 
Table H. This shows that the dynamics of the many- 
body system is also chaotic in the region of phase space 
that corresponds to the collective motion. Table || dis- 
plays a hierarchy of well separated rates A > 7 > a with 
A a. This hierarchy is particularly pronounced for 
large N since A (a) increases (decreases) with increas- 
ing N. The damping and equilibration rates differ for 
the following reason. While the dipole mode is damped 
out once the single-particle momentum distribution ex- 
hibits spherical symmetry, equilibration requires also a 
considerable change in the radial momentum distribu- 
tion. Note that the dipole moment D and the moments 
I v of the single-particle energy distribution effectively 
test the angular and radial momentum distribution, re- 
spectively. Our results show that a generic self-bound 
many-body system with two-body interactions displays 
a fast equilibration in the angular components opposed 
to a slow equilibration in the radial components of the 
single-particle momenta. This observation is also of in- 
terest with view on heavy ion collisions. 

To examine the role of the two-body interaction more 
closely, we consider a system of N independent parti- 
cles moving inside a circular billiard of diameter a. This 
system is integrable since single-particle energies and an- 
gular momenta are conserved quantities. It can be taken 
as a mean-field approximation of our Hamiltonian |2C| ] 
and is motivated by the observation that the surface of 
the interacting many body-system (i.e. the points in 
two-dimensional configuration space where interactions 
occur) becomes sharper with increasing numbers of par- 
ticles and may roughly be approximated by a circle of 
diameter a. As before, we take initial conditions at ran- 
dom from the region that corresponds to collective mo- 
tion and compute the dipole-moment (|J) as a function 
of time. Fig. ^ compares the evolution of the longitudi- 
nal dipole-moment D y obtained for the mean-field model 
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with the result obtained for the interacting many-body 
system. One observes damped oscillations with a period 
and a damping rate that are slightly larger than for the 
interacting system. The mean-field Hamiltonian there- 
fore captures the damping of collective motion quite ac- 
curately. In the integrable system, damping is due to the 
dephasing of single-particle oscillations that is induced 
by the collisions with the surface, i.e. the directions of 
single-particle momenta get randomized on a time scale 
I/7. The vadility of the mean- field picture is also con- 
firmed by the observation that the periods and damp- 
ing rates display no significant TV-dependence. This is, 
however, different for the Lyapunov exponents and equi- 
libration rates. These do depend on N through the pres- 
ence of the two-body interaction. Note that our observa- 
tion of N- independent frequencies to = 2tt/t and damp- 
ing rates 7 leads to a simple scaling law. Keeping the 
average single-particle energies E/N fixed and scaling 
the diameter of the iV-body system as a a N 1 / 3 yields 
ix>,7 oc TV -1 / 3 . These results are in semi-quantitative 
agreement with experimental data o, i.e. Tilo cx TV -1 / 3 
in heavy nuclei, and 7 oc N~^^ 3 '" 2 ^. 

We also considered initial conditions of the interact- 
ing many-body system with larger momentum spread 
Aq x w Aq y w q, i.e. the initial momentum distribu- 
tion is closer to an equilibrium distribution. This situ- 
ation is closer to the perturbative excitation of the nu- 
clear giant resonance. As expected, one finds a collective 
oscillation of smaller amplitude and comparable period, 
showing that the observed phenomena are robust. The 
determination of the damping rate is, however, difficult 
since amplitide and statistical fluctuations are roughly of 
the same magnitude. 

In summary, we have considered damping and equili- 
bration of collective motion in a self-bound iV-body sys- 
tems with two-body (surface) interactions. This system 
is predominantly chaotic and exhibits damped collective 
motion that leads to equilibration. There is a hierar- 
chy of three well separated time scales starting with the 
onset of classical chaos at short times, damping of the 
dipole mode at intermediate times and equilibration at 
long times, respectively. The damping is mainly due to 
dephasing and may be understood quite accurately in a 
mean-field picture of non-interacting particles. Conse- 
quently, periods and damping rates show no significant 
iV-dependence. Equilibration, however, requires the ran- 
domization of the magnitudes of single-particle momenta 
and a two-body interaction. Lyapunov exponents and 
equilibration rates depend mildly on N . The presented 
model exhibits a rather homogeneous phase space struc- 
ture and the main characteristics of self-bound many- 
body systems like atomic nuclei. Our results show that 
it captures important features of collective motion in such 
systems. 
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nucl- 



N 


3 


100 


150 


220 


330 


500 


A 


0.963 


1.506 


1.561 


1.62 


1.67 


1.75 


AA 


0.004 


0.017 


0.020 


0.02 


0.02 


0.02 


Runs 


70000 


100 


100 


50 


30 


10 



TABLE I. Lyapunov exponents (mean values A and RMS 
deviations AA) for iV-body systems with unit single-particle 
energy E/N = 1. Results are obtained from the listed number 
of runs, A is given in units of h/ma . 
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N 


100 


150 


220 


330 


500 


A 


1.52 ±0.02 


1.57 ±0.01 


1.61 ±0.03 


1.67 ±0.02 


1.73 ±0.02 


7 


0.39 ±0.05 


0.38 ±0.08 


0.42 ±0.06 


0.39 ±0.03 


0.37 ±0.04 


a 


0.21 


0.19 


0.15 


0.13 


0.10 


T 


1.77 ±0.05 


1.77 ±0.04 


1.76 ±0.03 


1.77 ±0.02 


1.77 ±0.01 



TABLE II. Damping rates 7, periods r, equilibration rates 
a and Lyapunov exponents A of the collective mode. (7,-7-, A: 
Mean values and RMS deviations are obtained from ten runs 
for each listed ./V-body system, a: exponent obtained from 
average of many runs.) The single-particle energy is E/N = 1; 
r, 1/A, I/7 and l/a are given in units of ma 2 /%. 
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FIG. 1. Time evolution of the dipole-moment (|3|) for a 
system of N = 220 nucleons at energy E — N. (Full 
line: longitudinal D H -component; dashed line: transversal 
Dx-component.) Time and dipole-moment are given in units 
of ma 2 jTt and a, respectively. 
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FIG. 2. Second moment !?,{£) of the single-particle energy 
compared to its equilibrium value 1% f° r N = 220. Time in 
units of ma 2 /%. 
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FIG. 4. Time evolution of the longitudinal dipole-moments 
D y for a system of N = 500 nucleons at energy E — N . 
(Full line: system with two-body interaction; dashed line: 
mean-field model.) Time and dipole-moment are given in 
units of ma 2 /h and a, respectively. 



12 3 

IPx.yl 




12 3 

|p| 

FIG. 3. Top: Normalized momentum distributions of single 
particle momenta \p x \ (top: full line), \p y \ (top: dotted line) 
and |p| (bottom: full line) after the equilibration of the col- 
lective motion compared to Maxwell-Boltzmann distributions 
with temperature kT = E/N = 1 (dashed lines). A system of 
N — 220 nucleons is considered; momenta are given in units 
of h/a. 
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